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Abstract: We introduce a new global pressure formulation for immiscible 
three-phase compressible flows in porous media which is fully equivalent to the 
original equations, unlike the one introduced in [S]. In this formulation, the total 
volumetric flow of the three fluids and the global pressure follow a classical Darcy 
law, which simplifies the resolution of the pressure equation. However, this 
global pressure formulation exists only for Total Differential (TD) three-phase 
data, which depend only on two functions of saturations and global pressure: 
the global capillary pressure and the global mobility. Hence we introduce a 
class of interpolation which constructs such TD-threc-phase data from any set 
of three two-phase data (for each pair of fluids) which satisfy a TD-compatibility 
condition. 
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Une formulation des ecoulements triphasiques 
en pression globale et totalement equivalente 

Resume : Nous introduisons unc nouvellc formulation en pression globale 
pour les ecoulements triphasiques immiscibles compressibles en milieu poreux. 
A la difference de [5] , cette formulation est totalement equivalente aux equations 
physiques originales, mais elle permet toujours d'exprimer le flux volumetrique 
total des trois fluides a I'aide de la pression globale par une loi de Darcy clas- 
sique, ce qui simplifie la resolution de I'equation en pression. Cette formula- 
tion n'existe cependant que pour des permeabilites triphasiques satisfaisant une 
condition de Differentielle Totale, qui ne dependent plus que de deux fonctions 
des saturations et de la pression globale : la pression capillaire globale et la 
mobilite globale. C'est pourquoi nous introduisons une famille d'interpolation 
permettant de construire de telles permeabilites a partir de la connaissance de 
trois jeu de donnces diphasiques (un pour chaque paire de fluides) vcrifiant une 
condition de compatibilitc. 

Mots-cles : Ecoulement triphasique en milieu poreux. 
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1 Introduction 

The numerical simulation of three-phase immiscible compressible flows in porous 
media requires the knowledge of three-phase relative permeabilities and capillary 
pressures. In practice, experimental values for these data are available only 
on the boundary of the ternary diagram, i. e. for the three two phase flows 
corresponding to each pair of fluids. The three-phase relative permeabilities are 
then derived by interpolation from these sets of two-phase data (see for example 
[5]). The existence of various interpolation formula is in itself a sign that none 
of these formula is the ultimate truth. 

We take advantage in this paper of this lack of experimental three-phase 
data, and introduce a new class of TD -interpolations, which is designed to 
simplify the numerical simulation of the flow by allowing the use of a global 
pressure formulation. 

The global pressure reformulation of the original flow equation was intro- 
duced for incompressible two-phase flows in [3l[T], and generalized, under the 
approximation that the volume factors be evaluated at the new global pressure 
instead of the corresponding phase pressure, to compressible two and three-phase 
flows in [S]. For three-phase flows (both incompressible and compressible), the 
original global pressure formulation required that the thec-phase relative per- 
meabilities and capillary pressures satisfy a Total Differential (TD)- condition. 
An algorithm for the determination of such TD-three-phase data was given in [S] 
and implemented in [H [7] , but its complexity has limited the use of the global 
pressure in numerical simulation codes. Nevertheless, comparison with other 
approaches [S] show the computational effectiveness of this approach when it 
can be put to work, which may explain the current revival of interest for the 
global pressure : the approximation on the volume factors required in the orig- 
inal global pressure formulation of compressible flows has been lifted recently 
for two phase flows in P] in an indepcndantly work, and is lifted for three-phase 
flows in the present paper. 

So we first show that, under an integral Total Differential (TD)- condition 
linking three-phase relative permeabilities and capillary pressures, the total vol- 
umetric flow of the three phases is governed by a single phase Darcy-like law 
for a new global pressure variable P. This new formulation is fully equivalent to 
the original three-phase compressible equations, in opposition to the formula- 
tion proposed in . Three-phase relative permeabilities and capillary pressures 
which satisfy the TD-condition we be called TD-three-phase data. 

Next we show that any collection of three two-phase data sets (defined on 
the boundary dT of the ternary diagram T) which is the trace of some TD-three- 
phase data (defined on T) satisfies necessarily a TD-compatibility condition. 

Finally, we study the TD-interpolation of a collection of three two-phase 
data sets which satisfies the TD-compatibility condition : we show that TD- 
interpolation reduces to the choice of two functions over the ternary diagram 
(one global capillary pressure, to be defined in scction|31 and one global mobility) 
which satisfy boundary conditions determined by the three given two-phase data 
sets. This is expected to make more easy the numerical determination of TD- 
three-phase data. 
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2 The three-phase immiscible compressible equa- 
tions 

Let the fluids be numbered in order of decreasing wettability, and denote by 
upper case letters Sj ,Pj,... saturation and pressure distributions (function of 
the space and time variables x,t)^ and by lower case letters Sj,pj, . . . saturation 
and pressure levels (real positive numbers). It will be convenient to use vector 
notations for the saturations : 



(1) 



(2) 



S = {81,83) = vector of saturation fields , 
s = (sijSa) = vector of saturation levels , 

and to denote by T and dT the ternary diagram and its boundary : 

T = {s such that < < 1 , j = 1, 3 , si + S3 < 1} 
dT = {s G T such that si = or S3 = or si + S3 = 1} 

2.1 Conservation laws : 

For each phase j = 1 (water), 2 (oil), 3 (gas) one has : 

^{(/)(x,Ppore)i3,(P,)5,}+V-^, =0 , J = 1,2,3. (3) 

where: 

X = vector of spatial coordinates, 

4>{x,p) = porosity at location X and pressure p, (4) 
Ppore = pore pressure, ~Pi,P2, or P3. 

and where, for each phase j — 1,2,3 : 

Pj ~ Pj{x,t) = pressure, 

8j = 8j{x,t) = reduced saturation , 8 = 81,82,83 

0<8j <1 , 81+82 + 83 = 1 , 
(fj = ifj{x,t) = volumetric flow vector at reference pressure. 

-^jiPj) ~ PjiPj)/pT^ = volume factor at pressure Pj, 
where : PjiPj) ~ density, {mass per unit volume) at pressure pj 

= density at reference pressure, 

(5) 

2.2 Muskat law : 

The volumetric flow vector of each phase j = 1, 2, 3 at reference pressure is given 
by: 

= -K{x) d,{P,) kr,{S){VP, - Pj{Pj)gVZ) (6) 

where: 

K{x) = absolute permeability at location x, 
'ijiPj) ~ Bj / fij = phase mobility at pressure Pj, 
f^jiPj) ~ phase viscosity at pressure Pj, 
krj{s) = phase relative permeability at saturation levels s, 
g = gravity constant , 

Z{x) = depth of location x. 
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2.3 Capillary pressure law : 



P1-P2 = P^^iSi) 
P3-P2 = Pf(53) 



(8) 



where: 



Pc'^i^i) — water — oil capillary pressure at water saturation level Si, 
Pc^i^s) — gas — oil capillary pressure at gas saturation level S3. 

(9) 



3 Classical resolution : "pressure equation" 

The numerical resolution of equations (I3|)([6|)(l8|) is usually done by solving one 
"pressure equation" , to be defined in the next section, with respect of one of 
the phase pressure, say P2 , and two of the "saturation equations" ^ , say with 
respect to Si and 5*3. 



3.1 Forming the "pressure" equation : 

Summing up the three conservation laws ^ gives : 

d ^ 

-{0(X, Ppore) J2 B^^Pl) 5,} + V • g = , (1) 

where q is the global volumetric flow vector: 

3 

9 =' E = -K^{^P2 + h^Pc^ + h^Pc' ' pg^z} . (2) 

i=i 

Here A, /i,/3,p are the global mobility, the water and oil fractional flows and 
the global density expressed as function of the oil pressure level p2, using the 
capillary pressure laws ([5]) : 

X{s,p2) = kri{s) di{p2 + P^\si)) + kr2{s) d2{p2) + krsis) d3{p2 + P!\s3)()3) 
f,{s,p2) = kr,{s)d,{p2+Pp{s,))/Xis,p2) , , (4) 

h{s,P2) = 1 - ,fl{s,P2) - f3is,P2) , (5) 
Pis,P2) = E /i(«'P2)/3jfe) • (6) 

In this approach, equations ([ijlll]) have to be solved for the oil pressure P2 at 
each time step.. 



3.2 Difficulties : 

Equation ((T|)([5|) is not a classical pressure equation, because of the gradient of 
capillary pressure terms in ([2]), which makes its discretization delicate. More- 
over, the individual phase pressures can be singular, as it is shown on figure [U 
for a two-phase problem : the difference between water and oil pressure is the 
capillary pressure, whose derivative is infinite at residual saturations - hence the 
water and oil pressure cannot be both regular across a water-oil front. So the 
choice of P2 as numerical unknown is likely to require a fine or adaptive mesh. 
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3.3 Let us have a dream . . . 

Could it be possible to find a "global pressure" field P such that : 

• P is smooth, in opposition to the individual phase pressures, 

• P governs the global volumetric flow of the three fluids : 

q = -Kd{VP-pgVZ} , (7) 

where now d and p denote the global mobility and density at a global 
pressure level p - compare with equation Q © : 

dis,p)'^= X{s,p2) , p{s,p)'^= p{s,p2) , (8) 

(we have kept for simplicity the same notation for the global density as function 
of p2 and p) 

4 Searching for a global pressure P 

Comparison of ([2]) ([7]) suggests to search for a new dependant global pressure 
variable P by setting : 

P^P2+Pcg{S,P) (1) 
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where the global capillary pressure function : 

{s,p) ^ F,gis,p) (2) 

is to be chosen such that : 

• ... P is in the range of the phase pressures. With the chosen wet- 
tabihty conventions, the capillary pressures ^ satisfy : 

PP{si)<0 , PP{1)^0 , ^(si)>0 



dsi 

dS3 



(3) 



so that 



-Pmin < -Pi < -P2 < -Pj < Pmax at all poiuts X and time t , (4) 

where Pmin is a known lower bound to the water pressure Fi , and Pmax is 
a known upper bound to the gas pressure P3. 

One wants that the global pressure satisfies a similar inequality : 

Pmin < Pi < P < P3 < Pmax at all poiiits X and time t . (5) 

This will make P a natural candidate for the pore pressure Ppore required 
to evaluate the porosity </> (see 

«... our dream comes (almost . . . ) true ! Comparison of equations 
(HI and ([7]) suggests to require that : 

for any saturation/global pressure fields 5*1(0;, t), S3{x,t), F{x,t), 
the global capillary pressure function ^ satisfies: 

VF,g{S, F) = MS, F - F,g{S, P)) VP,i2(5i) + his, F - F,g{S, P)) VPf (Sg) 



P2 P2 
+ ^(5,P)VP, 

(6) 

where V denotes, as in equations (H])© the gradient with respect 
to the space variables x. When this condition is satisfied, equation ^ 
becomes, with the notations ^ : 

f)P 

q = -Kd{{l - ^)VP - pg\/Z} , (7) 

which is the desired equation ([7|) up to the factor (1 — dFcg/dF) ! 

• ... P is uniquely determined by ([T]) once P2 and S = [Si, 6*3) are 
known. This property will make it possible to use (P, S) as dependant 
variables instead of (P2 , S) . It will be satisfied as soon as the stability 
condition : 

QP 

\-^{s,p)\ < 1 for all S e T , P^in <P< Pmax (8) 
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is satisfied. Notice that the computation of P from S, P2 is required only 
for the determination of the initial global pressure distribution; after that, 
during the course of the computation, the numerical code knows only P, 
and it is P2 which is computed from P and S using ([1]). 

Hence a global pressure formulation will exist as soon as one can find a global 
capillary pressure function s,p Pcg{s,p) which satisfies the three conditions 



4.1 Satisfying condition (Ej) : 

Condition ^ will constrain the global capillary pressure up to a constant, so 
we can fix its value at one point of the ternary diagram. Choosing : 

Pcg{l, 0,p) for P,„i„ <p< P„,ax , (9) 

will imply the desired property ^ as it follows immediately from (|19p below. 



4.2 Satisfying condition (I6D : 

Taking successively Si{x,t) = constant and S3{x,t) = constant shows that ^ 
will be satisfied if and only if the global capillary function satisfies : 

dP dP^'^ 

-Tr^{s,P) = fl{s,p- Pcg{s,p))-T^{si) , 

kg, . „ , ..d% (10) 



(s,p) = h{s,p- Pcg(,S,p))—^{S3) , 



.S3 ds. 

for all S e T , P,„in <p< Pmax- 

This relation implies that, if it exists, Peg satisfies a differential equation 
along any smooth curve C of T. Hence the existence of a global capillary pressure 
function Peg satisfying (fT(I|) amounts to the following Total Differential (TD) 
condition : 

for any s G T, Pniin < P < ^'max and any smooth curve C : [0, 1] T, such 

that C(0) = (1,0) , C(l) = s, the solution of the differential equation : 

d3 dP'^'^ dP^"^ 

^ = /i(C,p-/3)^(Ci)Ci + /3(C,p-/3)^(C3)C^ , /3(0) = 
at as I as3 

satisfies (3{1) = s indepcndantly of C. 

(11) 

When this TD-condition is satisfied, the global capillary pressure satisfying 
© © is given by : 

Peg{s,p)=(3il)- (12) 



An alternative formulation of the TD-condition is obtained by requiring that 
the second derivative d^Pcg/dSidSs takes the same value when computed from 
either line of III 



There exists a function s,p Pcg{s,p) s.t. for all s G T , Pmin < P < -Pmax 

(13) 
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Figure 2: The Total Differential condition for three phase data : Pcg{s,p) is 
required to be independant of the curve C along which it is computed. 



In the approximate global pressure formulation (see [S]), the Peg terms were 
neglected in the right-hand sides of (fTU]) . and hence in So condition (fT^ 

could be used to check wether a given set of three-phase data satisfied the 
TD-condition or not, and to construct fractional flow functions satisfying this 
condition [4]. This approach cannot be used anymore in the exact formulation 
given here, where the condition (|13|1 begins with "there exists Peg such that 
..." : one has to resort to condition pip to check if a given set of three-phase 
data satisfy the TD-condition, and to use a more direct - and hopefully simpler 
- approach for the construction of TD-three-phasc data, which we dccribe now. 

The TD-condition pip or (fT3|) has the effect of reducing the number of 
functions which can be chosen freely over the ternary diagram: instead three 
relative permeabilities kri , kr2 , kr^ which are function of saturations (or equiv- 
alently two fractional flows /i , /a and the global mobility A which are function 
of saturations s = (si, S3) and oil pressure P2), one can choose only in a global 
pressure formulation the global capillary pressure Peg and the global mobility d 
which arc function of saturation s ~ (si, S3) and global pressure p. 

Based on ^TU\i the fractional flows and relative permeabilities associated to 
Peg and d are : 

uj{s,p) = dPeg/ds,{s,p)/ dPl^/ds,{sj) , J = 1,3 (14) 
V2{s,p) = l-vi{s,p)-v^{s,p) (15) 
fcr,(s,p) = p,{s,p)d{s,p)/dj{p~Peg{s,p) + Pr-{s,)) J = 1,3, (16) 
kr2{s,p) = {1 - iyi{s,p) - 1^3(8, p))d{s,p)/d2{p - Peg {s,p)) , (17) 

where we have used the notation i/j for the fractional flows considered as function 
of the global pressure p. They are related to the fractional flows fj as function 
of the oil pressure p2 by: 

Vo{s,p) = fj{s,p- Peg{s,p)) , j = 1,2,3. (18) 
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4.3 Satisfying condition 

Let now the three-phase data satisfy the TD-condition pT|) . Pcg{s,t) be the 
global capillary pressure defined by ([5]) (fTT|) (fT^ . s = (51,53) £ T, P,„i,i < p < 
Pmax be given, and : t {1 — t + tsi,0) and : i ^ (si, ts^) be the curves 
shown in figure [31 The functions P^{t) = Pcg{C^ {t),p) , j = 1,3 satisfy the 
differential equations (sec (fTTj) ) : 



12 



-E- = -f,il-t + tsu0,p-(3^)^{l-t + ts,){l~s,) , /3i(0) 

^ = +/3(5i,ts3,p-/33)--^(ts3)s3 , /33(0) = /3i(l) . 
at as3 



(19) 




Si water si oil 

Figure 3: One convenient path for the computation of dPcg/dp. 

Derivation of with respect to p shows that the functions 7^ {t) — dPcg/dp{Cj (t) , p) 
are given by the the differential equations : 

Wpl2 

-il~t + tsi,0,p- /3I) —^(1 - t + isi) (1 - si)(l - 7I) , 7^(0) = , 
asi 

= +^(^1: i^3,P - P') ^(tss) ssil - 7') , 7^(0) = 7^(1) . 

ds3 



' ^7! 




dt ~ 


dp2 


d-1^ _ 




dt ~ 


dp2 



Both equations in (|20p are of the form 

dj 



dt 
so that 



a(i)(l-7) , 7(0) = 70 < 1 



(20) 

7(f) = l-(l-7o)exp{- / a(r)dr} , 

(21) 



BP 

^{s,p) = 7'(1) = 1 - ^M-J^ {<^\r) + a\T))dT} 



(22) 



which shows that the coefficient l~dPcg/dp in the global Darcy law ^ is always 
strictly positive. 
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The stability condition (0) that \dPcg/dp\ < 1 is then equivalent to : 

Vs e T , Vp G [Pmin, Pmax] , exp{- / (a^r) + a3(T))dT} < 2 . (23) 

Jo 

Derivation of the fractional flows /i(s,p2) and /3(s,p2) given by ©(HI) with 
respect to p2 gives the following expressions for and (remember : kr^ ~ 
on Ci) : 

d' d' 

kridikr2d2{-j- - ^) ,pi2 

^ ^ {kridi + kr2d2Y dsi ^ ' ' ^ ' 

d' d' d' d' 

krsdskndii-^ - -J-) + kr3d3kr2d2{-^ - -j-) ,p32 

a3(t) - + ^3 rfi 'ds d2' dF^ 

''^*) - + (^kr,d,+kr2d2 + kr,d,Y ds, '^^^> 

where krj,dj,d'j , j = 1, 2, 3 and , j = 1, 3 are evaluated at saturations and 
phase pressures associated to the point C^{t) (for ai) or C'^{t) (for 02) of T at 
the given global pressure level p. 
When the fluid densities satisfy : 



max J 



Vs G T , VpG [Pn.in,P, 

fip- P,,is,p) + Pf (.3)) >f(P- P^^^'^P)) ^ T^P P^^^'^P) + ^c'(^i)) 

«3 «2 "1 

(26) 

a^{T) and Q;^(r) are positive, so that : 

dP 

VS G T , Vp G [Pmin, Pmax] , < -g^ (s , p) < 1 , (27) 

and the stability condition ^ is satisfied. 

Naturally, is only a sufhcicnt conditions for In practice, implemen- 
tation of the global pressure formulation requires the determination of dPcg/dp 
over the ternary diagram T and the global pressure range of interest [Pmin, Pmax], 
for example by formula (P^(|24p(P5|) . so one can check directly on d Peg /dp if 
condition ([S]) is satisfied. In the case where condition @ is not satisfied, it is 
not known wether this limitation is due to the global pressure formulation, or 
is the sign that the original equations are not anymore well-posed, or change of 
type. 

Remark 1 In the case of a water-oil-gas system, one has d'l ~ d'2 ~ 0, and 
condition i26]) is satisfied : the global capillary pressure, when it exists, satisfies 
always ([27]) • In that case, dPcg/dp = on the water - oil side ofT, so that 
= and h22\l simplifies somewhat. H 
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5 TD-interpolation of two-phase data given on 
the sides of the ternary diagram T 

Let two-phase relative permeabilities and capillary pressure curves be given on 
the boundary dT of T for each system of two-fluids : 

water ~ oil : si ^ fcrP(si), A:r2^(si) P^^{si) , 

water - gas : sj fcrP(si), krl^{si) Pc^{si) = Pc^(si) - ^^^(1 - si) , 
gas -oil : S3 ^ krl^iss), krj^iss) P^^iss) , 

(1) 

as well as the three phase mobility functions (see ([T])) : 

P,-^d,{pj), J = 1,2,3, (2) 

The question arises of wether one can find, for each given global pressure 
level Pmin l£ P ^ -Pmax, thrce-phase relative permeabilities s G T ^ krj{s,p) G 
M, j = 1, 2, 3 (compare with (O) which . . . 

1. honor the two-phase data on OT, 

2. and satisfiy the TD-condition PT|) on T. 

Satisfying condition 1 and 2 amounts to search for a global capillary pressure 
function {s,p) Pcg{s,p) and a global mobility function (s,p) d{s,p) such 
that the fractional flow Vj{s,p) defined by pO|) : 

dP,g . . I dP^^ , , 

osi I ds\ /ON 

dP^g. , / dPf. , 

produce three-phase relative permeabilities krj{s,p) : 

''^^^^'^^ = rfr(p-pji,pKP-(.r)) 

d2(p-^'cg(s,p)) 

'^'^'^'^ ^ d3ip-p.2s^p)\ppis.)) '^'^'^^ 

which coincide with the given two phase data ^ on dT. 

5.1 Determination of Pcg{s,p) on dT. 
Let Pmin < p < Pmax bc givcn, and denote by 

C^^(i) = (l-t,0) , C^^{t)^{0,t) , C^^(t) = (l-t,t) , 0<t<l(5) 

a parameterization of the three edges of dT. When it exists, the global capillary 
pressure [t) = Pcg{C^^ {t)),i,j — 1,2,3, z < j on dT satisfies, according to 

m ■■ 
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• along the water-oil-gas sides 



dt 
dt 



= - - i, 0, p - (1 _ t) , (0) 



dsi 



dcf 12,datE 
= V. 



jp32 



dS3 



dcf 23. data/ 
= Vi 



.3 '(t,p) 

• along the water-gas side : 



/r-(i-t,i,p-/?^^) 



12 



■> dsi 



il-t) 



clef 13, data 



jp32 



^ ds?, 



dot 13.dat£ 



I /3"(0) 



, 



(6) 



(7) 



where ff^'^'^{s,p2) and f^^^^{s,p2) denote, for s e OT, the two-phase fractional 
flows as function of the oil pressure derived from the data ((T|) ^ using two phase 
versions of ©([I]). 

Using the TD-condition (jlip . one sees that the existence of a global cap- 
illary function implies that the two-phase data ([I]) ([2]) satisfy necessarily the 
TD-compatibility condition 



(8) 



i.e., using ©([T]) and the two-phase fractional flows 
defined there : 



12. data 23. data 13, data 13, data 
)i'3 ''^l '1^3 



12, data 13, data' 



dP. 



12 



dsi 



23, data 13, data 



for P,„i„ <p<P„ 
(9) 

This compatibility condition is a constraint on weighted means of the frac- 
tional flows, and hence a non-linear constraint on the mean values of the two- 
phase relative permeability data sets ([T]). When it is satisfied, one sees that the 
global capillary function Peg has to satisfy the Dirichlet boundary condition : 



P — P' 



.data d?if 
eg 



(water-oil edge), 
(water-gas edge), 
/3^^ (gas-oil edge). 



(10) 



5.2 Determination of 



eg , 



dn 



s,p) on dT. 



Now that Peg is known on dT, equation (fT(I|) gives, for s G dT 



dsi 
dPcc 



dP^^ 

/f-*-(s,p-P,t*"(s,p))-^(si) 



ft'%S.P~Pcg'%S,p)) 



ds 
dPl 



ds:f 



-is,) 
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Let then n be the unit outer normal to T. The normal derivative of Peg on dT 



is given by : 



dn 
dn 
dn 



_ n/3 



- 2- 



9s9 



■) (water-oil edge) 



\/3, 



dsi 



V3(dPcg 



dsi 



) (water-gas edge) , 
-) (gas-oil edge) , 



(12) 



Combining pT|) and (fT^ shows that Peg satisfies the Neumann boundary 
dition: 



dPcg 

dn 



ap 



data 



dcf 



^3 ^12, data 



dP, 



12 



9n 



V3 



13. data 



3 "1 
dP}2 



IS.data'^^c^N 



dsi 

V3 23, data "J^e 
dS3 



dS3 



(water-oil edge), 
(water-gas edge), 
(gas oil edge). 



1 12, data 23. data 13. data 13, data 

where ' , , , ly^ ' 

Conversely, let Peg satisfy the Dirichlet and Neumann conditions ((TU]) (dH). 
Then, on the water-oil edge for example, (fTU]) and the first line in ^ implies 



are defined in ([51) 0. 



(13) 



that ui{C'^'^{t),p) = V- 



12. data 



(t,p), and US]) nil) and the first fine of ((T3l) implies 



that j/3(C^^(t),p) = for < i < 1. Hence the fractional flows associated to a 
global capillary function which satisfies (fTUl) honor the data on dT. 



5.3 Determination of d{s,p) on dT. 

Now that the global capillary pressure P^^*^ is known on dT, honoring the three 
sets of two-phase relative permeabilities ((l} on dT amounts simply to impose 
on d the Dirichlet boundary condition : 



where : 



jdata 



on dT , 



(14) 



jdata 



dcf 



kri di{p 
fcrP di{p - 
fcr|^ ^3 {p - 



P^^'^ + PP) + kr^^ d2{p 

pdata I pl2\ 



eg 



+ krl'^ dsip 

.23 (p. 



kr^-' d2 



pdata'^ 
eg ^ 
pdata 1 

CQ ' 

^data'^ 



eg 



P?) 



(water-oil), 
(gas- water) , 
(gas-ofl). 

(15) 



5.4 Determination of Peg and d on the interior of T. 

Conditions PH)) (|13p p3)) are the sole conditions to be satisfied to ensure that the 
TD-threc-phase data derived from Peg and d match the three given sets (HJ® 
of two-phase data on dT - provided the latter satisfy the TD-compatibility 
condition ([9]). 

There is hence a large choice of functions Peg and d to choose from to perform 
TD-interpolation. 
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For example, one can use the smoothest functions which satisfy the boundary 
conditions dTO]) (US]) (HH) by setting: 



These equation can be solved by finite element over the ternary diagram T. 

One can also try to match some a-priori given three phase target permeability 
model on the interior of T, knowing that the match will not be exact, unless the 
target relative permeabilities happen to satisfy the TD-compatibility condition 
pip . This requires the implementation of an optimization algorithm. 

6 Conclusion 

1. When the three-phase relative permeabilities and capillary pressure satisfy 
the TD-condition ill]) , a global capillary pressure function Pcg{s,p) exists 
such that the classical compressible immiscible three-phase flow equations 
are fully equivalent to a global pressure formulation : the total volumetric 
flow q of the three phases and the global pressure P = P2 + Pcg{S,P) 
follow the global Darcy law ([7]) . 

2. In the global Darcy law ([7|), one has always 1 — dPcg/dp{s,p) > 

3. If the compressibility of the fluids increases when their wettability de- 
creases (condition (dH)), one has always 1 > 1 — dPcg/dp{s,p) > 0, and 
the stability condition 0) is satisfied, which ensures that one can deter- 
mine P from P2 and S. Condition ([26]) ) is always satisfied in water-oil-gas 
flows where the compressibility of water and oil is neglected. 

4. TD-three-phase relative permeabilities can be obtained by interpolation 
of three sets of two-phase data on the sides of the ternary diagram, pro- 
vided the two-phase data satisfy the TD-compatibility condition 0). TD- 
intcrpolation amounts to choose a global capillary function Peg and a total 
mobility function d which satisfy boundary conditions (fTU|) (fO|) but 
can be chosen freely inside the ternary diagram. This interpolation class 
takes advantage of the lack of information on actual three-phase data in- 
side the ternary diagram to simplify the numerical simulation of the flow. 
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